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Since the 1970’s with the establishment of the big bang model, it has become clear that 
some of the most restrictive constraints on certain neutrino properties come from astro- 
physical and cosmological considerations. Furthermore, in 1987 the detection of neutrinos 
from the supernova in the Large Magellanic Cloud provided a new “neutrino laboratory” 
as well as confirming our basic understanding of gravitational collapse energetics. This ar- 
ticle will review those constraints on neutrinos derived from cosmological and astrophysical 
considerations. 

We will first review the freeze out of neutrinos in the early Universe and derive the 
cosmological limits on masses for stable neutrinos. We will then use the freeze out argu- 
ments coupled with observational limits to constrain decaying neutrinos as well. We will 
also review the limits to neutrino properties which follow from SN19S7A. We will then 
look at the constraint from big bang nucleosynthesis on the number of neutrino flavours. 
Before ending, we will briefly look at astrophysical constraints on neutrino- mixing as well 
as future astronomical observations of relevance to neutrino physics. 

Cosmological Mass and Decay Limits 

Cosmological limits to neutrino mass and decay properties depend on their relic number 
density from the early Universe. If a massive particle species remained in thermal equilib- 
rium until the present, its abundance, n/s ~ (m/T) 3 ^ 2 exp(— m/T), would be absolutely 
negligible because of the exponential factor (s = entropy density). If the interactions of the 
species freeze out (i.e., T < H) where T is the interaction rate and H is the cosmological 
expansion rate at a temperature such that m/T .is not much greater than 1, the species 
can have a significant relic abundance today. W T e will now calculate that relic abundance. 

First, suppose that the species is stable (or very long-lived compared to the age of the 
Universe when its interactions freeze out). Later we will consider the case where the species 
is unstable. Given that it is stable, only annihilation and inverse annihilation processes, 
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e-g., 

vv * — >XX, (1) 

can change the number of n’s and P’s in a comoving volume. Here A' generically denotes 
all the species into which n’s can annihilate. In addition, we assume that there is no 
asymmetry between v's and P’s. 

We will also assume that all the species A', X into which v, 9 annihilate have thermal 
distributions with zero chemical potential. Because these particles will usually have addi- 
tional interactions which are ‘stronger’ than their interactions with n's, the assumption of 
equilibrium for the A”s is almost always a good one. For example, let A’, X = e~, e + ; 
while the neutrinos only have weak interactions, the e^s have weak and electromagnetic 
interactions. 

The evolution of the number density can be expressed 1 in terms of the total anni- 
hilation cross section (cr 4 |r' | ) 

~^f + 3 Hn » = -{^h’Dh^-C 7 ^) 2 ] ( 2 ) 

This equation for the evolution of the abundance of a species is a particular form of 
the Riccati equation, for which there are no general, closed-form solutions. Before we 
solve the equation by approximate methods, let’s consider the qualitative behavior of the 
solution. The annihilation rate varies as n^Q times the thermally-averaged annihilation 
cross section (cr^ |) . In the relativistic regime, (rn^/T <C 3) heq ~ T 3 , and like other 
rates, will vary as some power of T. In the non-relativistic regime, {m v /T ;§> 3) 
n EQ ~ (mT) 3 / 2 exp(— m/T), so that decreases exponentially. In either regime, 1P,4 
decreases as T decreases, and so eventually annihilations become impotent , roughly when 
T A ~ H , which we call freeze out. 

• Hot Relics: First consider the case of a particle species such that rn/T < 3 at freeze 
out. In this case, freeze out occurs when the species is still relativistic and the equilibrium 
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number density per comoving volume Y E q = ;? E q/s is not changing with time. Since > E q 
is constant, the final value of I' (Y = nfs) is very insensitive to the details of freeze out, 
and the asymptotic value of V, Y(m — » oo) = 1^,, is just the equilibrium value at freeze 
out : 

3oc — ^eq = 0.27S 5 eff/^5. (3) 

where g t ^ — g (bosons), 0.75g (fermions), and g counts the internal degrees of freedom. 
Thus the species freezes out with order unity abundance relative to entropy $ (or the num- 
ber density of photons). Assuming the expansion remains isentropic thereafter (constant 
entropy per comoving volume), the abundance of i/s today is (s 0 is the present entropy 
density) 


= 2970F OC cm 3 

( 4 ) 

825 {g^/g.) cm -3 . 

( 5 ) 


If, after freeze out, the entropy per comoving volume of the Universe should increase, say 
by a factor of 7, the present abundance of z/s in a comoving volume would be diminished 
by 7. 

A species which decouples when it is relativistic is often called a hot relic. The present 
relic mass density contributed by a hot relic is simple to compute: 

p„ — soYr^m — 2.97 x lO^o^m/eV) eV cm -3 , (6) 

Q u h 2 = 7.83 x KT^eff/^.s] (m/eV). (7) 

Based upon the present age of the Universe we know that Qoh 2 < 1; applying this bound 
to the contribution of the species v to Oo h 2 we obtain a cosmological bound to the mass 
of the u: 

m < 12.SeV[<7„s(.T / )/ g t n] (S) 
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Light (mass < MeV) neutrinos decouple when T ~ few MeV, and g.$ = g, = 10.75. 
For a single, 2-component, neutrino species g e R = 2 x (3/4) = 1.5, so that g e n/g.s = 0.140. 


This implies that 


nji 2 

m„ 


__ m v 
91.5 eV’ 
< 91.5 eV. 


( 9 ) 

( 10 ) 


This cosmological bound to the mass of a stable, light neutrino is often referred to as the 
Cowsik-McClelland bound. 2 (In their original paper, Coswik and McClelland consider a 
4-component neutrino (g = 4), and took 12 < 3.S, h = 1/2 and T„ = T, which resulted in 
the bound m < 8 eV.) 

If there are more than one light (< MeV) species, this bound applies to the sum of the 
masses of the light neutrinos. 


• Cold Relics: Now consider the more difficult case where freeze out occurs when the 
species is non-relativistic ( m/T ^ 3), and I E q is decreasing exponentially with m/T. In 
this case the precise details of freeze out are important. 

First we will parameterize the temperature dependence of the annihilation cross sect ion. 
On general theoretical grounds we expect the annihilation cross section to have the velocity 
dependence oc v 1 ’ . where p = 0 corresponds to 5-wave annihilation, p = 2 to p - wave 

annihilation, etc. Since (v) ~ T 1/2 , (<7^|r|) oc T n , n = 0 for s-wave annihilation, n = 1 for 
p - wave annihilation, etc. Therefore we parameterize (o^M) as 

(a A \v\)=o 0 (T/m) n (11) 

With this parameterization, the Boltzmann equation for the abundance of r/’s becomes, 

dY/dx = -Xx- n ~ 2 (Y 2 - l' E 2 Q ), (12) 


where 


x = m/T. 


(13) 
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(14) 


A = 0.264(g.s/gl /2 )mp L ma o , 

1’eq = 0.145(g/g mS )x 3 ' 2 c~ r . (15) 


As shown in ref. 1, eq. 12 can be solved approximately to good accuracy where it is 
found that 


_ 3.79(7} + l).r" +1 

(g m s/gl^)m PL mcr 0 

where 


(16) 


x f ~ /??[0.03S(?7 + 1 )(g/gl /2 )m PL mcr 0 ] (17) 

-(?7 + 1/2)/?? /??[0.03S(?? + l)(<//t/y 2 ) 777p£,?77<7 0 ] 

As with a hot relic, the present number density and mass density of cold relic j/s is 
easy to compute, 


??i /0 = s 0 Ioo = 29701’^, cm 3 

(77 + l)x} +1 


1.13 x 1(T 


{y*s/gl /2 )m PL mo 0 
^ ( ?7 + l)x” +1 Ge\ 


ft./? 2 = 1.07x10 


cm 


-3 


7-1 


( 18 ) 

(19) 


[g*s / 9 l * 2 ) m PLV 0 

(where the subscript / denotes the freeze out value). It is very interesting to note that the 

relic density of ids is inversely proportional to the annihilation cross section and mass of 
the particle 


1' = 


3.79(77 + l){gV 2 /g.s)x f 


mm PL {a A \v\) ' ^ 

The smaller its annihilation cross section, the greater its relic abundance — the weak prevail. 

Moreover, the present mass density only depends upon the annihilation cross section at 
freeze out, which for 7 ? = 0 (s-wave annihilation) is independent of temperature (and 
energy). 
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Let us now look at the specific application of this to massive neutrinos (n? MeY). 

Annihilation for such a species proceeds through Z° exchange to final states il; where 
7 — v Li e, r > u,d, s, ... ( v L denotes any lighter neutrino species). The annihilation 
cross section depends upon whether the heavy neutrino is a Dirac or Majorana species; for 
T < m < Mz, the annihilation cross section is 


Mv\) 


Dirac 


G 2 F m 2 

2tt 


£0 -zf^Kc^ + cl) 


+ -Z?(C 2 , + C%)} 


foM)l 


Majorana — 




+cl 2z 2 ] 
^ .4 1 ? J’ 


( 21 ) 


where = m-Jiv, ,3 is the relative velocity, and C\ and C'.\ are given in terms of the 
weak isospin J3. the electric charge q. and the Weinberg angle 0\y by C , j = y'3. Cy = 
J’A — 2</sin 2 6\y. (Ye have assumed that the neutrino is less massive than Mz-) 

In the Dirac case, annihilations proceed through the 5-wave and {<7.4 |r|) is velocity 
independent: 

o 0 ~ c 2 G 7 F m 2 /2-K (22) 

where c 2 ~ 5. Taking <7 = 2 and g, ~ 60, from our formulae we find 


xj ~ 15 + 31n(m/GeV) + ln(c 2 /5) 


r “ “ 6 * 10 ' 8 G^ 


-3 


31n(m/GeV) ln(c 2 /5) 
1 d ~ h 


15 


15 


(23) 


from which we compute that 


n^h 2 = 3 (m/GeV)" 2 

where we have included the identical relic abundance of the antineutrino species (Cl„p = 
20„). Note that freeze out takes place at T) ~ m / 15 ~ 70 MeV(m/GeY) - before the 
interactions of light neutrinos freeze out. This is because as neutrinos annihilate and 


1 + 


3ln(m/GeVj 


15 


(24) 
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m(eV) 

Fig. 1: The contribution to fi 0 ^ 2 for a stable neutrino species of mass m (from ref. 1). 

become rare, the annihilation process quenches. Requiring Q „/? 2 < 1 we obtain the so- 
called Lee- Weinberg bound: 

m > 2 GeV. (25) 

Although it is often called the Lee-Weinberg bound, the basic argument 3 was noted a 
decade earlier by Zeldovich, Novikov and Chiu. 

For the Majorana case, annihilation proceeds through both the s and p- waves; however 
the formulae which obtain for xj , and f } l/ h 2 are similar. In Fig. 1 we show the con- 
tribution to Q. 0 h 2 for a stable, massive neutrino species. For m < MeV, Q t ,h 2 ocmas the 
rehc abtmdance is constant. For m > MeV, f l^h 2 oc m~ 2 as the relic abundance decreases 
as m -3 . The relic mass density achieves its maximum for m ~ MeV. 

Neutrino masses less than about 92 h 2 eV, or more than about 2 GeV (Dirac) or about 
5 GeV (Majorana) are cosmologically acceptable. 

These limits are quite impressive when compared with the laboratory limits, u t , at 
250 keV and v 7 at 35 MeV, and imply that both must be below 92 eV if they are stable. 
Furthermore, recent searches for the products of neutrino annihilations in the sun and earth 
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by Kamiokande and Irvine-Michigan-Brookhaven (IMB) probably constrain any stable 
massive neutrino to be < 12GeV or the high energy neutrino^ produced by annihilations 
would have been observed. 4 

Before leaving stable massive neutrinos, it is worth noting that they can still be the 
dominant mass in the Universe. Relic neutrinos of a few GeV mass provide closure den- 
sity and behave as cold dark matter. Moreover, this possibility may soon be tested by 
more sensitive searches for their annihilation products, and/or cryogenic detectors. Relic 
neutrinos of mass ~ 30 eV provide closure density and behave as hot dark matter. While 
laboratory experiments will eventually probe a v t mass as small as 10 eV, we will probably 
have to wait for the next nearby supernova to probe iq* and v r masses in the 30 eV range. 
While hot dark matter and adiabatic density perturbations (such as those produced by 
inflation) seem to be incompatible with observations, hot dark matter with cosmic strings 
(as the seed perturbations) is a very viable and interesting structure formation scenario. 

Unstable u's 

Now consider the possibility of an unstable neutrino 5 species whose decay products are 
relativistic, even at the present epoch. It is clear that the mass density bound for such 
a species must be less stringent: from the epoch at which they decay (say, ~ = Z£>) until 
the present, the mass density of the relativistic neutrino decay products decreases as J? -4 , 
as opposed to the i? -3 had the neutrinos not decayed. Roughly speaking then, the mass 
density today of the decay products is a factor of (1 + X£>) -1 less than that of a stable 
neutrino species. 

The precise abundance of the neutrino decay products is very easy to compute. Denote 
the energy density of the relativistic decay products by po, and for simplicity we will 
assume that they do not thermalize. The equations governing the evolution of the daughter 
products axe: 

Pd + 4 Hp D = pjr. 
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(26) 


P»( R ) = Pu{ R ,) (j^j exp(-</r). 

where R t , t, is some convenient epoch prior to decay, t, < r. The relic density of the 
decay products is obtained by integrating (26): 

PD(t) = PviT - 1 jf ~R^ exp(— t'/r )(H'. (27) 


Assuming that around the time the neutrinos decay (t ~ r) the scale factor R oc i n 
(n = 1/2 radiation dominated; n = 2/3 matter dominated) we can evaluate this integral 
directly, and find that the present density of relic, relativistic particles from neutrino decays 
is 

H(t) 

PD(to) = n\p„{to)—=r^ (28) 

where p^(to) is the present density that neutrinos and antineutrinos would have had they 
not decayed, and R(r) is the value of the scale factor at the time t — t. As expected, 
the present energy density of the decay products is less than that of a stable neutrino 
species, by a factor of n!J t ?(r)/i? 0 ~ (1 + zd) 1 - During the matter-dominated epoch 
(t > 4.4 x 10 lo (fi 0 fc 2 ) -2 sec), R(t)/R 0 = 2.9 x 10 _12 (Q o ^ 2 ) 1 ^ 2 /c 3 i so that the reduction 
factor is 

n\R(r )/R q = 2.6 x 10 _12 (D o /i 2 ) 1/3 r s 2 / c 3 . (29) 

During the radiation-dominated epoch, R(t)/Ro = 2Ag7 1/12 x 10 _10 tj/ 2 , so that the reduc- 
tion factor is 

n!i?(r)//2o = 2.1 x 10' 10 9 ; 1 / 12 t^ c 2 . (30) 

Using the results of our earlier calculations for we obtain the following constraint 
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-4 -2 0 2 4 6 0 10 12 

log[m(eV)] 

Fig. 2: The forbidden region of the neutrino mass-lifetime plane based upon the requirement that 
Q„h 2 < 1 (from ref. 1). 

to the epoch of decay (for neutrino masses which fall in the previously disallowed range) 

m< 4 x 10 n eV(/! /12 r“ c ,/2 (light, r < < EQ ) 

m < 4 X 10 53 eV(fi o /j 2 )- 1/3 r- c 2/3 (light, r > t E q) 

m > 3 x 10- 5 GeV( ? 7 1/2 V s 1 e / 4 (heavy Dirac, r < t EQ ) 
m > 3 X 10~ 6 GeV(Q o h 2 ) 1 / 6 T*/ c 3 (heavy Dirac, r > # EQ ) ' ; 

m > 7 X 10 _5 GeA r ^r 1//24 t s (/ c 4 (heavy Majorana, r < * E q) 
m > 8 x 10 -6 GeV(Q o /) 2 ) 1/,6 r s ] e { 3 (heavy Majorana, r > t E g). 

The excluded region of the neutrino mass-lifetime plane is shown in Fig. 2. (Considera- 
tion of the formation of structure in the Universe leads to a significantly more stringent 
constraint to the mass density of the relativistic decay products; structure cannot grow in 
a radiation-dominated Universe. For a discussion of these constraints see ref. 6.) 

The limits just discussed 5 apply irrespective of the nature of the decay products (so 
long as they are relativistic). If the decay products include “visible” particles, e.g., photons, 
e ± pairs, pions, etc., much more stringent limits can be obtained'. We will now consider 
the additional constraints which apply when the decay products include a photon. (For 
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the most part these same limits also apply if the decay products include e* pairs.) The 
limits which obtain depend both qualitatively and quantitatively upon the decay epoch, 
and we will consider five distinct epochs. 

Before discussing these limits, it is useful to calculate the time at which the energy 
density of the massive neutrino species would dominate the energy density in photons. 
The energy density in photons is = (7r 2 /15)r 4 , and assuming the neutrinos are NR, 
their energy density is p „ = Taking g,s — 4, the energy densities are equal when 

T ~ 3 Y^.m. For heavy neutrinos Y. x is given by eq. 20, and for light neutrinos. Y x , ~ 0.04. 
Thus we find that the relic neutrino energy density will exceed the photon energy density 
at T/m < 0.1 for light neutrinos, and T/m < 2 x 10 -8 m^ v for heavy neutrinos. Using 
t — 1 sec/Tj^y for the age of the Universe, the epoch of matter domination (by massive 
neutrinos) is given by 


f(sec) 


light neutrinos 
heavy neutrinos. 


(32) 


10 14 (?7}/leY) 2 
3 X 10 9 77?Q eV 

(Here, and throughout the following discussion, “light' 1 will refer to neutrinos of mass less 
than an MeV, and “heavy” will refer to neutrinos of mass greater than an MeV, but less 
than Mz-) 


• tu — 3 x lO^sec < r: If the neutrino lifetime is greater than the age of the Universe, 
neutrinos will still be decaying at the present and decay-produced photons will contribute 
to the diffuse photon background. Assuming that the neutrinos are unclustered (the most 
conservative assumption), the differential number flux of decay -produced photons (per cm 2 
sr sec erg) is 


cLF y ri v c 

dEdQ 4 tttH 0 


E~ l 


3/2 


, m / 2 , 


(E < m/2) 


(33) 


where for simplicity we have assumed that each decay produces one photon of energy 
m/2 and that Q, 0 = 1. Taking the number flux to be dT^/dYl ~ EdE 1 /dEdQ. and H 0 = 


12 



log(E/eV) 


_ e _ 6 -4 -2 0 2 4 6 8 10 12 14 16 18 20 



Fig. 3: The diffuse photon background. Vertical arrows indicate upper limits, and horizontal arrows 
indicate integrated fluxes (> E.) Circles and triangles indicate the total cosmic ray flux (p's, nuclei, and 
photons) which provides an absolute upper limit to the photon flux at. the highest energies (from ref. 1). 

50km sec -1 Mpc _1 , we find 


d 

dn 


~ 10 29 r sec cm 2 sr 1 sec 1 light neutrinos 
~ 3 x IO^t^ttjqpV cm -2 sr -1 sec -1 heavy neutrinos 


(34) 


A summary of the observations of the diffuse photon background are shown in Fig. 3. 
The differential energy flux, dT / dEdEl, is shown as a function of energy and wavelength. 
From this data, a very rough limit of 


d!F. y f MeV ^ _ 2 _i _j 

dn~ l^J cm sr sec ’ 


(35) 


can be placed to the contribution of neutrino decay-produced photons to the photon back- 
ground. Based upon this, the following lifetime limit results: 


Tsec > 


10 23 m e v light neutrinos 
10 25 m5ev heavy neutrinos, 


(36) 


applicable for neutrino lifetimes r •> 3 x lO^sec. The forbidden region of the mass-lifetime 
plane is shown in Fig. 4. 
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• t rec ~ 6 x 10 ]2 (ft 0 ft 2 r 1/2 sec < t < t v : If neutrinos decay after recombination, but 
before the present epoch, then the decay-produced photons will not interact and should 
appear today in the diffuse photon background. Again, for simplicity, assume that each 
neutrino decay produces one photon of energy m/2. Then the present flux of such photons 
is 

dT^ n„c 

dQ 47r 

~ 3 x 10 11 cm -2 sr -1 sec -1 light neutrinos 
— 4 x 10 4 rn^ev cm -2 sr -1 sec -1 heavy neutrinos (37) 

Where we have assumed that when the neutrino species decays, it is non-relativistic, so 
that each decay-produced photon today has energy E ~ m/2(l -f ~o)~ where (1 + : D ) ~ 
3.5 x 10" (n 0 h 2 r 1/3 r-c /3 - Comparing these flux estimates to our rough estimate of the 
diffuse background flux we obtain the constraints, 

m < 2 x 10 6 (f lo/? 2 ) -1 / 3 ^// 3 eY light neutrinos 

m > SxlO- 3 (^ 2 ) 1/f V s ^ 3 GeV heavy neutrinos, (38) 

applicable for neutrino lifetimes in the range 3.5 x 10 11 (f2 o ^ 2 ) _1 ^ 3 sec < r < 3 x 10 17 sec. For 
very light neutrino species the assumption that the species decays when it is non-relativistic 
breaks down. If the species decays after t = and before the present epoch, and is 

relativistic when it decays, the decay-produced photons will be comparable in energy and 
in number to the CMBR photons and will cause significant distortions to the CMBR. Thus 
a neutrino species which decays while relativistic in the time interval 10 6 ^<3x 10 17 
sec is forbidden. The excluded region is 200 < t sec /m eV <4x 10 2O (f2 0 /i 2 ) 1//3 , for 


f 3.5 x 10 8 (ft o fc 2 )- 1/3 f- 2 / 3 t stc > 4.4 x 10 lo (fU 2 r 2 
eV ~ \ 4.6 x 10 H-W t sec < 4.4 x 10 l0 (Qoh 2 )- 2 

The forbidden region of the mass-lifetime plane is shown in Fig. 4. 


(39) 


• Uherm — 10 6 sec < t < t T(C : For neutrino decays which occur during this epoch, the 


14 




Fig. 4: Cosmological limits to the mass and lifetime of an unstable neutrino species which decays radia- 


tively (from ref. 1). 
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decay-produced photons can scatter with electrons, which can in turn scatter with Cosmic 
Microwave Background Radiation (CMBR) photons, thereby changing the spectral shape 
of the CMBR 7 . However, during this epoch processes which can alter the number of 
photons in the CMBR, e.g., the double Compton process, 7 + e — ► 7+7 + e, are not effective 
(i.e., T < H). Therefore, the result of dumping significant amounts of electromagnetic 
energy density from neutrino decays is a Bose-Einstein spectrum (with p-, ^ 0) for the 
CMBR. The CMBR is to a very good precision a black body. Thus, any electromagnetic 
energy density resulting from neutrino decays during this epoch must be much less than 
that in the CMBR itself. Recalling that 

p„ _ mlocS 

Pi Pi 

~ 0.1 m/T light neutrinos 

~ 2 x 10 -s m5evWT heavy neutrinos, (40) 

and requiring that p ( ,/p-y ^ 1, we obtain the following limits for a neutrino species that, 
decays during this epoch: 

777 < 10' r- c 1/2 eY light neutrinos 

777 > 4 x 10“ 3 GeV heavy neutrinos, (41) 

where we have taken t sec — . These limits are applicable for neutrino lifetimes in the 

range 10 6 sec < r < 10 13 sec. The forbidden region of the mass-lifetime plane is shown in 
Fig. 4. (A neutrino species which decays after nucleosynthesis and produces photons of 
energy greater than 30 MeV can lead to photofission of the light elements produced during 
nucleosynthesis; additional, more stringent bounds result 9 ). 

• tend n udeo — 3 min < t < t tfie rm- For neutrino decays which occur during this epoch, 
the decay- produced photons can be thermalized into the CMBR because both Compton 
and double Compton scattering are effective (T > H). However, in so doing the entropy 
per comoving volume is increased. This has the effect of decreasing the present value of r t 
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relative to the standard scenario. It is known that luminous matter (necessarily baryons) 
provides Q.lum ~ 0.01, and thus provides direct evidence that today ij>4x 10 -n . On 
the other hand, primordial nucleosynthesis indicates that at the time of nucleosynthesis 77 
corresponded to a present value of (3 — 10) x 10~ 10 (ref. 8). Thus any entropy production 
after the epoch of nucleosynthesis must be less than a factor of ~ 10 _9 /4 x 10 -11 30 . 

This leads to the limits 


10 9 

> 

m eV uVc 2 

light neutrinos 


10 7 

> 

m G>V Tig 

heavy neutrinos, 

( 42 ) 


applicable for neutrino lifetimes in the range 200 sec < r < 10 6 sec. This bound too is 
shown in Fig. 4 (also see refs. 10). 

• tbtgin nucho — lsec < t < t tnd nuc i (0 : If the neutrino lifetime is longer than about a sec. 
then massive neutrinos can contribute significantly to the mass density of the Universe 
during nucleosynthesis, potentially leading to an increase in 4 He production. Only the 
equivalent of 1 additional neutrino species can be tolerated without overproducing 4 He. 
One additional neutrino species is about equivalent to the energy density contributed by 
photons. Since the crucial epoch is when the neutron-to-proton ratio freezes out (/ ~ 1 
sec, T ~ 1 MeV), the constraint that follows is (p v / Pi)r~Me\ ^ 1 - This results in the mass 
limit 

m > 5 x 10 _3 GeV heavy neutrinos. (43) 

Note there is no corresponding limit for a light species because a light species is just one 
additional relativistic neutrino species. This limit, which is applicable to a heavy neutrino 
species with lifetime greater than about 1 sec is shown in Fig. 4. 

• t« 1 sec: A neutrino species which decays earlier than about 1 sec after the bang 
disappears without leaving much of a cosmological trace. Its decay products thermalize 
before primordial nucleosynthesis, and its only effect is to increase the entropy per comov- 
ing volume. If we understood the origin of the baryon-to-entropy ratio in great detail, and 
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could predict its ‘pre-nucleosynthesis’ value, then we could use entropy production by the 
decaying neutrino species to obtain constraints for very short lifetimes. 


• Astrophysical implications: Neutrino decay into visible modes can have “astrophys- 
ical” effects too (refs. 11). As the detection of neutrinos from SN 19S7A dramatically 
demonstrated, type II supernovae are a copious source of neutrinos. The integrated flux of 
neutrino-decay-produced photons from type II supernovae that have occurred throughout 
the history of the Universe can be used to obtain a very stringent bound to acceptable 
neutrino masses and lifetimes. 


Each type II supernova releases about 3 x 10 53 ergs of energy in thermal neutrinos 
with average energy about 12 MeV — or about A T „p ~ 5 x 10 5 ' neutrinos and antineutrinos 
of each species. The historical (last 1000 yrs) type II rate in our own galaxy is about 
1 per 30 yrs (give or take a factor of 3), and the observed extragalactic rate is roughly 
l.l/? 2 per 100 yrs per 10 10 Lb®. Using the measured mean blue luminosity density of the 
Universe, Lbq ~ 2.4/? x 10 8 £b,t, Mpc -3 , this translates into a present type II rate (per 
volume) of Tsn — 2.5 /? 3 x 10~ 85 cm -3 sec -1 . Assuming that the type II rate has been 
constant over the history of the Universe (a bold assumption), the differential photon 
number flux is 




9 TsNtlN^ 


(44) 


d9.dE 5v / 2 47r(E (/ )r/m (E„y/*EW 

where for simplicity we have assumed that the supernovae neutrinos are mono-energetic, 

with = {£„) ~ 12 MeV, that each decay-produced photon carries half the energy of the 
parent neutrino, and a flat Universe. For this energy spectrum (£-,) = (E^/6 ~ 2 MeV. 
Comparing the expected photon number flux at 

1 T SN tlN uP m 
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with the measured diffuse 7-ray flux at a few MeV, 3 x 10 -3 cm -2 sr _I s -1 , we obtain the 
following constraint: 


fsec >5x 10 12 (IW3 x 10 8o cm 3 sec 1 )??? e v- 


(46) 
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Of course, this bound only applies to neutrino species light enough to have been produced 
in supernovae (m < 10 MeV) and which decay outside the envelope of the exploding star 
( r sec 10 -5 m e v) by the present epoch ( t < 10 n m e v sec). (Based upon 7-ray observations 
of SN 1987A made by the SMM spacecraft a similar, slightly more restrictive bound ob- 
tains. Furthermore, the lack of observed ionization around SN 19S7A by v — » v t + e“ limits 
this mode for m > 1 MeV (see refs. 11).) This constraint is shown in Fig. 5. 

For a neutrino species which decays within the envelope of the exploding star, and 
thereby deposits energy in the envelope a different bound can be derived. Any energy 
deposited by neutrino decays in the envelope will be thermalized and radiated in the 
visible part of the spectrum. The energy radiated by SN 1987A in the visible was only 
about 10 4 ' ergs, while each neutrino species carries off about 10 53 ergs! The energy which 
is deposited in the envelope by a hypothetical, unstable neutrino species is 

Edep - min[l, R B sg/~lab ] (47) 

— mi n [ 1 0 53 ergs , 1 0 48 7?7 eV / r sec ergs] (48) 

where Rbsg ~ 3 x 10 12 cm is the radius of the envelope of the progenitor blue super giant 
(Sanduleak -69 202, by name), and r LAB = {E^r/tn is the neutrino lifetime in the rest 
frame of the supernova. Comparing this to the observed energy of 10 4 ‘ ergs, we obtain the 
bound 


™ e vMec <0.1 

(t aec > 10 5 m eV ) 

(49) 

m eV > 10 T 

(t 3ec < 10 -5 m eV ) 

(50) 


This constraint too is shown in Fig. 5. 


A neutrino species which can decay radiatively, iq — > v, + 7, necessarily has an elec- 
tromagnetic coupling that may be quantified as a transition magnetic moment, //,j = 
K,j(e/2m e ). The transition magnetic moment and neutrino mass and lifetime are related 


19 


6 


3 

0 

-3 



-4 -2 0 2 4 6 8 10 12 

log[m(eV)] 


Fig. 5: Astrophysical limits to the mass and lifetime of an unstable neutrino which decays radiatively 
(from ref. 1). 





by 


t 1 = QEMK 2 m 3 /8m 2 e , 

k = 0.44 T~] /2 m~y /2 , (51) 

where we have assumed 777^ m t . The transition moment leads to an electromagnetic 
correction to v — e scattering. Laboratory limits to v — e scattering through the transition 
moment leads to the bound < 10" 8 , or 

r sec > 2 x 10 15 77 ?e - x 3 (i/ /( Vt + 7). (52) 

Further, such a transition moment leads to neutrino pair emission from white dwarfs and 
red giants through the process plasmon — + v t v v For k ~ lO -10 — 10 -n plasmon vv emission 
can be a very significant cooling mechanism for these objects, and can effect their evolution. 
Based upon this, a limit of H tJ < 10“ 10 or so has been derived for neutrinos less mass than 
10 keV (see, e.g., the paper of Beg, et al. 1 1 ) . This translates to the limit 

r sec > 2 x 10 19 t777v (™ £ 10 keV). (53) 

All of the astrophysical and cosmological constraints just discussed are summarized in 
Figs. 4 and 5. These constraints serve to illustrate how a large variety of cosmological and 
astrophysical observations can be used to probe particle properties in regimes beyond the 
reach of the terrestrial laboratory. 

Limits To the Number of Families 

Another area where cosmological and astrophysics constraints have been important is in 
limiting the number of neutrino families, A 7 ^. These arguments take on added importance 
now that accelerator experiments are beginning to check them 12 . The most important 
bound comes from big bang nucleosynthesis 13 . A second, very different but less stringent 
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bound, comes from SN19S7A H , Let us first look at the big bang nucleosynthesis argument. 

The power of big bang nucleosynthesis comes from the fact that essentially all of the 
input physics is well determined. The relevant temperatures, 0.1 to 0.05 MeV, are well 
explored in nuclear physics labs. Thus, what nuclei do under such conditions is not a 
matter of guesswork but is precisely known. In fact, the nuclear physics is known far 
better for these temperatures than it is known in the centers of stars like our sun. The 
temperature at the center of the sun is only a little over 1 keV. This energy is below the 
energy where nuclear reaction rates yield significant results in laboratory experiments, and 
only the long times and higher densities available in stars enable anything to take place at 
all! Unfortunately, for stellar astrophysics this means that nuclear reaction rates must be 
extrapolated to many orders of magnitude below their laboratory-measured values. The 
big bang laboratory does not have this problem. The reactions occur at temperatures and 
densities where cross sections and the like are known and well studied in the laboratory. 

To calculate what happens, all one has to do is follow a gas of baryons with density 
pb as the Universe expands and cools. As far as nuclear reactions are concerned, the 
important epoch begins a little above 1 MeV and ends a little below 100 KeV. At higher 
temperatures, no complex nuclei other than single neutrons and protons can exist, and 
the ratio of neutrons to protons, n/p, is just determined by thermodynamic equilibrium, 
n/p = e~Q ! T , where Q = 1.3 MeV is neutron-proton mass difference. Equilibrium applies 
because the weak interaction rates are much faster than the expansion of the Universe at 
temperatures much above 1 MeV. At temperatures much below 0.1 MeV, the electrostatic 
repulsion of nuclei prevents nuclear reactions from proceeding as fast as the cosmological 
expansion separates the particles. 

After the weak interaction drops out of equilibrium, around 1 MeV, the ratio of neu- 
trons to protons changes more slowly, by free neutrons decaying to protons and similar 
transformations of neutrons to protons via interactions with the ambient leptons. By the 
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time the Universe reaches 0.1 MeV, the ratio is slightly below 1/7. For temperatures above 
0.1 MeV, the high entropy of the Universe suppresses the abundance of nuclei. Once the 
temperature drops to about 0.1 MeV, nuclei begin to be present in significant amounts, 
starting with 2 D adding neutrons and protons, making 3 H and 3 He. These, in turn, capture 
neutrons and protons to produce 4 He or 3 H and 3 He can collide to also yield 4 He. Since 4 He 
is the most tightly bound nucleus (in this region of the periodic table), the flow of reactions 
converts almost all the neutrons that exist at 0.1 MeV into 4 He (for neutron/proton ratios 
less than unity). The two-body chain essentially ceases there, because there are no stable 
nuclei at either mass-5 or mass-S. Since the baryon density at big bang nucleosynthesis is 
relatively low (much less than 1 g/cm 3 ) only reactions involving two-particle collisions 
occur. It can be seen that combining the most abundant nuclei neutrons, protons, and 4 He 
via 2-body interactions always lead to unstable mass-5. Even when one combines 4 He with 
rarer nuclei like 3 H or 3 He, we still only get to mass-7, which when hit by a proton, the 
most abundant nucleus around, yields mass-S. Eventually, 3 H radioactively decays to 3 He, 
and any mass-7 made, radioactively decays to ‘Li. Thus, big bang nucleosynthesis makes 
4 He with traces of 2 D, 3 He, and 'Li. (Also, all the protons left, over that did not capture 
neutrons remain as hydrogen.) All other chemical elements are made later in stars and in 
related processes. (Stars jump the mass-5 and -S instability by having gravity compress 
the matter to sufficient densities that 3-body collisions can occur and jump the mass-5 
and -8 gaps.) A neutron/proton ration of ~ 1/7 yields a resultant 4 He primordial mass 
fraction, Y = 2(n/p)/(n/p+ 1) m 0.25. 

The only cosmological parameter in such calculations is the density of the baryon gas at 
a given temperature. From the thermodynamics of the expanding Universe we know that 
Pb oc T 3 , thus we can relate the baryon density at 10 U K to the baryon density today, when 
the temperature is about 2.75 K. The problem is, we don’t know pi, today, so the calculation 
must be carried out for a range in pb- The cosmological expansion rate depends on the 
total mass-energy density. For cosmological temperatures much above 1 eV the energy 
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density of radiation exceeds the mass-energy density of the baryon gas. Thus, during big 
bang nucleosynthesis, we need the radiation density as well as the baryon density. The 
baryon density determines the density of the nuclei and thus their interaction rates, and the 
radiation density controls the expansion rate of the Universe at those times. The density 
of radiation is just proportional to the number of “types'' of radiation. Thus, the density of 
radiation is not a free parameter provided we know how many types of relativistic particles 
exist at temperatures ~ 0.1 - 1.0 MeV. 

Assuming that the relativistic particles at 1 MeV are photons, c. //, and t neutrinos 
(and their antiparticles) and electrons (and positrons), the big bang nucleosynthetic yields 
have been calculated for a range in present pb (more precisely the baryon to photon ratio), 
going from less than that observed in galaxies to greater than that allowed by the observed 
large-scale dynamics of the Universe. The 4 He yield is almost independent of the baryon 
density, with a very slight rise in the density due to the decreasing entropy per baryon, 
which enables nucleosynthesis to start slightly earlier, when the neutron/proton ratio was 
higher. No matter what assumptions one makes about the baryon density, it is clear that 
4 He is predicted by big bang nucleosynthesis to have to be around 25% of the mass of the 
Universe, This was first noted by Hoyle and Tayler 15 and later found by Peebles 16 and by 
Wagoner, Fowler, and Hoyle 1 '. The current results do not differ in any qualitative way 
from Wagoner, Fowler, and Hoyle’s original detailed calculations. 

The fact that the observed helium abundance in all objects is about 20 to 30% is 
certainly a nice confirmation of these ideas. Since stars produce only a yield of 2% in all 
the heavy elements combined, stars cannot easily duplicate such a large 4 He yield. While 
the predicted big bang yields of the other light elements were also calculated in the 1960's, 
they were not considered important at that time, since it was assumed in the 1960's that 
these nuclei were made in more significant amounts in stars. 18 However, work by our group 
at Chicago 8 , and others, thoroughly established big bang nucleosynthesis and turned it 
into a tool for probing the Universe, by showing that other light element abundances had 
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major contributions from the big bang and that the effects of any stellar contributions 
could be removed by appropriate techniques. Today the big bang predictions for all four 
light isotopes are used to test the model and use it as a probe of conditions at early times. 

In particular, it was demonstrated in the early 1970's that contrary to the ideas of 
the 1960's, deuterium could not be made in any significant amount by any realistic con- 
temporary astrophysical process 19 . The big bang deuterium yield decreases rapidly with 
increasing p b . At high densities deuterium gets more completely converted to ‘‘He; quanti- 
tatively this means that the present density of baryons must be below ~ 5 x 10" 31 g/cm 3 in 
order for the big bang to have produced enough deuterium to explain the observed abun- 
dance. Similar, though more complex, arguments 20 were also developed for 3 He, and most 
recently for 'Li, so that it can be said that only if the baryon density is between 2 x 10~ 31 
g/cm 3 and 5 x 10 -31 g/cm 3 are all the observed light element abundances consistent with 
the big bang yields. If the baryon density were outside of this narrow range, a significant 
disagreement between the big bang predictions and the observed abundances would result. 
To put this in perspective, it should be noted that for this range in densities, the predicted 
abundances for the four separate species cover a range from 25% to one part in ~ 10 10 . 
The big bang yields all agree with only one freely adjustable parameter, p< : . 

Recently, several non-standard scenarios of primordial nucleosynthesis have been proposed 
however, these scenarios with their additional adjustable parameters seem to be unable 
to account for the abundances of the 4 light isotopes, especially 7 Li. This speaks to the 
remarkable success of the standard scenario of big bang nucleosynthesis. 

This narrow range in baryon density for which concordance occurs is very interesting. 
Let us convert it into units of the critical cosmological density for the allowed range of 
Hubble expansion rates. From big bang nucleosynthesis 19 ’ 20 , it follows that the baryon 
density fig is less than 0.12 and greater than 0.03 (once one includes 21 age constraints on 
a flat Universe); that is the Universe cannot be closed with ba.ryov.ic ma.tt.er. If the Universe 



is truly at. critical density, then non-baryonic matter is required. This argument has led to 
one of the major areas of research at the particle-cosmology interface, namely, the search 
for non-baryonic dark matter. 

Another important conclusion regarding the allowed range in baryon density is that it 
is in very good agreement with the density implied from the dynamics of galaxies, including 
their dark halos. An early version of this argument, using only deuterium, was described 
over 15 years ago 22 . As time has gone on, the argument has strengthened and the fact 
remains that galactic dynamics and nucleosynthesis both suggest densities of about 10% 
of the critical density. Thus, if the Universe is indeed at critical density, as many believe, 
it requires that the bulk of the matter not be associated with galaxies and their halos, as 
well as being non-baryonic. 

With the growing success of big bang nucleosynthesis, the predictions came under more 
scrutiny. In particular, the 4 He yield was examined in detail since it is the most abundant 
of the nuclei, and thus in principle it is the one which observers should be able to measure 
to highest accuracy. In addition, it is very sensitive to the n/p ratio. 

In the standard calculation it is assumed that photons, electrons, and the three known 
neutrino species (and their antiparticles) are present in the Universe at the time of nucle- 
osynthesis. However, by doing the calculation with additional species of neutrinos we can 
see when 4 He yields exceed observational limits. The bound on 4 He comes from observa- 
tions of helium in many different objects in the Universe. However, since 4 He is not only 
produced in the big bang but in stars as well, it is important to estimate what part of the 
helium in some astronomical object is primordial, from the big bang, and what part is due 
to stellar production after the big bang. To do this we 24 have found that the carbon content 
of the object is well suited for tracking the additional helium produced. Carbon is made 
in the same mass stars that also produce 4 He, thus as the carbon abundance increases, so 
must the helium. (Other heavy elements such, as oxygen have been used previously, but 
these elements are not produced in the same mass stars as those that produce the bulk of 
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the helium.) The extrapolation of helium to zero carbon content in an object should be a 
good estimate of the primordial helium. We obtain ~ 0.235 as our best estimate for the 
mass fraction of helium produced in the big bang. The upper bound is what is important 
here. We formally estimate a three standard deviation bound as 0.247. In particular, it 
seems clear that the primordial 4 He was at least a little less than 25%. Since objects have 
heavy elements and possibly some associated extra-stellar produced helium and still have 
helium abundances of 25%, this certainly seems like a very safe upper bound. In fact, if 
anything our estimates are on the high side due to possible systematic errors, e.g., Pagel 25 
finds collisional excitation reduces the 0.235 to 0.233. 

We find (see Figure 6) that three (or two) types of neutrinos fit the data well, and a 
fourth is only marginally allowed if helium slightly exceeds the 3-a upper bound; any more 
neutrinos are strictly prohibited. Since each family contains a neutrino, we are saying that 
the total number of families is three or at most four. Thus, all the fundamental families of 
elementary particles may have been already discovered. Of course, this assumes that the 
neutrinos are “light,” i.e., less massive than ~ MeV. 

Supernova 19S7A and Neutrino Counting 

Let us now compare this bound with the supernova constraint. As is now well appreci- 
ated, neutrinos were detected 26 from SN 1987A by both Kamiokande 26 and IMB 27 . Both 
of these H^O detectors are most sensitive to u e + p — ► n . + e + because of its larger cross 
section. 

If the P e flux is assumed to come from a Fermi-Dirac (F-D) distribution at temperature 
T and total u e energy, e Pe , both IMB and Kamiokande are simultaneously fit with T ~ 4 
to 4.5 MeV and e Pc ~ 3 to 4.5 x 10 52 ergs. These figures are in remarkable agreement with 
the standard model 28 for gravitational core collapse of a massive star, if = 3. Thus, 
we have confidence that we have witnessed such a core collapse, and that we have a good 
understanding of its physics. Let us now turn the argument around and see how sensitive 
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our expected fluxes are to A r „. 

In a collapse to a neutron star, the binding energy, cb, must be radiated as neutrinos. 
The initial neutronization burst of i/ e ’s carries away a fraction /„ < 10% of eg on a timescale 
of < 10 ms. The remaining energy comes out in thermal vv pairs from reactions like 

e + e - — + vv (54) 


where through neutral currents all species of neutrinos with m„ < 10 MeV are emitted. 


Since electron scattering rates are small compared to z" e capture, even with five times 
more free electrons than protons, at most we expect one or two scattering events in the 
detectors for a SN at 50 Kpc (distance to LMC). Thus, the detectable fraction of is e Ce 
where 


2-Y„ 




(55) 


assuming an equipartition of energy emitted in the various neutrino species, as is found in 
the detailed models. (While average energy per neutrino is higher for and u T . their flux 
is correspondingly lower.) The number of counts, 7? , one expects in a detector of mass. 


M d , is 


(a) 2 Mb 


(56) 


(E Pt ) 4 t:R 2 18 m p 

where m v is the proton mass, R & 50 Kpc is the distance to LMC, {Ep c ) is the average 
P e energy, and (a) is the cross section appropriately averaged over a F-D distribution 
with appropriate threshold factors and efficiencies taken into account. The temperature of 
i/ e ’s is found to be ~ 3.2 MeV ( (£),) ~ 10 MeV) to good accuracy. Temperatures are very 
insensitive to model parameters being determined by microphvsics at the neutrinosphere 14 . 
The temperature for P e ’s is somewhat higher due to the smaller opacities at late times as 
protons disappear in the core, thereby minimizing charged current interactions. This 
enables the P e ’s to come from deeper in the star. Mayle et al. 28 find T Pc ~ 4 MeV in good 
agreement with the temperature inferred from the observations. (They do find a higher 
than thermal high energy tail to the distribution which can effect the high threshold IMB 
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but not Kamioka.) For detectors like Kamiokande where the threshold is well below the 
peak of the cross section weighted distribution, it is reasonable to use 


W) « Void + 'Mil H2T? 


(57) 


(For IMB a more careful procedure must be applied due to its high threshold.) Substituting 
into Eq(56) yields 
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which for Mp = 2.14 ktons (Kamioka) we obtain a prediction of 11 counts for N u = 3. 
While they actually observe 11, one should weigh their counts by efficiency effects to obtain 
16.5 ± 5. Solving for N„ yields 




(2 ± 0 . 6 ) 


\(JF>. t (_... ■> ) 

~ f n \ ( 50Kpcy 

. \ 4MeV J V2 x 10 53 ergs/ 

’V 0.9 ) V R )\ 


(59) 


Let us now see how high we can push this. While models can be found with f n > 0.1, 
it is obvious that 1 — /„ can never exceed unity. The effective T r/t . as used above, varies 
by < 25%. The binding energy for 1.4M S neutron stars (the mass of the collapsing core) 
is found to vary from 1.5 to 3 x 10 53 ergs for a wide range of equation-of-state 29 . Thus, 
we choose 3 x 10 53 ergs (4 xlO 53 ergs) as an (extreme) upper bound. The distance to the 
LMC varies in the astronomical literature by < 7%. We’ll adopt an extreme limit of 10% 
consistent with current SN 19S7A determination of the distance 30 . Combining all these 
extreme value yields 

< 6. 6(8. 9). (60) 


A more careful calculation taking into account different thresholds for both IMB and 
Kamiokande to obtain measured e Ce for predicted yields at. the T Pe inferred from the data 
yields essentially the same result (A T „ < 6. 7(9.0)) as given above. Thus, SN 19S7A gives 
a limit to comparable to accelerator experiments but not as strong as the big bang 
nucleosynthesis limits. 
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Fig. 6 : Helium mass fraction versus the baryon-tophoton ratio 77 . The lower bound of 2 x 10“ 10 derives 
from the 3 He 4- D and ( Li constraints, and the upper bound of 7 X 10 “ 10 from the D and 7 Li constraints. 
The three lines for each neutrino family correspond to neutron half-lives of 10.4, 10.5, and 10.6 minutes 
(from ref. 8 ). 




Other Constraints from SN19S7A 


SN19S7A has proven to be an amazing neutrino laboratory. In addition to the pre- 
viously mentioned limits, it has placed limits to the charge and magnetic moment of the 
neutrino that exceed current laboratory limits and its constraint to the mass of v t is com- 
parable to the best laboratory limits. Let us briefly review these bounds. 

Magnetic Moment 

Barbieri and Mohapatra, and Lattimer and Cooperstein 31 have shown that the obser- 
vation of P e ’s from SN 19S7A constrains the value of the magnetic moment of the neutrino 
to < 10 -11 //B- The argument is twofold, involving in a crucial way the fact the interaction 
cross section of right-handed neutrinos must be significantly weaker than those of left- 
handed neutrinos. (Right-handed Dirac neutrinos must interact more weakly so that they 
do not get counted in the big bang nucleosynthesis arguments 32 .) First, there is the limit 
from cooling the proto-neutron star too rapidly if z-7's can change to up's as a result of 
magnetic moment interactions in the proto-neutron star core. Second, there is the effect 
that a flipped up can escape from the higher temperature inner core and then get flipped 
back to a up by the intergalactic magnetic field. This latter situation could yield 70 MeY 
P £ ’s which were definitely not detected. It is argued that these processes limit the magnetic 
moment to < 10 -13 with 10" 11 //# as an extreme upper limit. However, Okun 33 has argued 
that these arguments can be circumvented if the magnetic moment is not static but is a 
majoron transition moment or if an appropriate MSW mixing 36 of neutrino species also 
occurs in the supernova. 


Neutrino Charge 


Barbelini and Cocconi 34 have argued that the absolute value of any neutrino charge 
must be < 10 -1 ' | e |; otherwise electrostatic repulsion would have spread the neutrino 
burst greater than ~ 10 sec on its 170,000 yr. flight from the LMC to earth. 
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Neutrino Mass 


Since the observed neutrino burst was relatively narrow (< 10 sec), despite energies 
which spanned a range of about a factor of two, it is obvious that any neutrino rest mass 
must be very small. While the relationship between mass, timespread and energy is a 
simple one, the key here is to decide on the significance of the time and energy spread, 
and to estimate what the intrinsic spread was in the neutrino burst in the absence of finite 
masses. 

The crucial, but simple, relationship at the heart of any analysis to constrain the u ( 
mass from the IMB and Kamiokande data is that for the time delay suffered by a neutrino 
during its flight to earth: 


. 1 R m 2 (m/10eV) 2 

At = - — — ~ 2.6sec — — — (61) 

2 c E 2 [Ej lOMeV) 2 S 

From this simple equation for At , it is clear that any mass constraint which follows will be 

in the general range of about 20eV, or so, which is comparable to existing laboratory limits. 

Given the sparseness of the data set (19 events in total), the subtleties of the detectors 

(response, thresholds, etc.) and the absence of a very specific, well-accepted standard 

model of the initial cooling, it is not surprising that many authors have ‘derived' limits 

(and even values!) for the mass ranging from a few eV to 30 eV. The most extensive 

and careful analyses to date 35 provide limits of around 20eV - 25 eV. While SN 1987 A has 

not really improved existing bounds, it is interesting that the constraint which is found is 

comparable to the present laboratory limits. 

Neutrino Mixing 


Neutrino mixing has been proposed as a solution to the solar neutrino problem 36 and 
the HomestaJke and Kamiokande observations of solar neutrinos place constraints on al- 
lowed mixing parameters 3 '. A supernova could potentially also test neutrino mixing 38 . If 
neutrino mixing occurs between supernova emission and detection, it can obviously alter 
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the detected neutrino signal. 

If MSW mixing is indeed the solution to the solar neutrino problem, then only u e 
mixing is possible and not v t -> v»(v T ). Thus, the solar neutrino solution would not 
affect the v e flux. However, it could deplete the initial neutronization burst. Unfortunately, 
there is no conclusive evidence that even a single v e + e~ — > u e -f c~ scattering event 
associated with the neutronization burst was seen. 

If we drop the solar neutrino solution and go to general MSW mixing, then we can mix 
v^{v T ) into P t , which might enhance the energy slightly, but would otherwise do little. No 
effect would occur for the electron scattering ’s. Thus no definitive statement can be 
made from SN 19S7A about neutrino mixing and oscillations. 

Secret Interactions 

Precious little is known about any interactions neutrinos may have beyond the standard 
weak interactions, e.g., additional neutrino-neutrino interactions as in the Majoron model. 
Since neutrinos from the supernova traversed 170,000 light years through the cosmic seas 
of relic neutrinos (and perhaps other particles such as majorons) without apparent atten- 
uation, any unknown (i.e., secret) interactions they might have with neutrinos (or other 
particles in the sea of relics) can be constrained: 

^secret ^ 10~ 25 CIU 2 . (62) 

Radiative Decays 

The fluence of neutrinos from SN1987A was enormous, ~ 10 10 cm -2 per species (in- 
tegrated over the observed burst). On the other hand there was no observation (above 
instrument background) of any high energy 7 rays: based upon the data of the Gamma 
Ray Spectrometer aboard the Solar Maximum Mission and 7 ray detectors on the Pio- 
neer Venus Orbiter a 7-ray fluence limit for the same time period of < 1 cm -2 follows. 
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This means that less than about 1 in 10 10 of the supernova neutrinos could have decayed 
producing a 7-ray. From these non-observations of 7-rays a limit of 

T »tc/ m t\' Z 2 x lO 15 ^ (63) 

can be set to the radiative decay of any neutrino species. Here B-, is the branching ratio 
for the radiative decay mode. 
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